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ABSTRACT 

Virtual effects of gravitons in the production of diphotons at the upgraded Tevatron and 
at the LHC are analysed with the idea of probing the parameter space of the Randall- 
Sundrum (RS) model. It is shown that this process yields stringent constraints on the 
parameter space of the RS model. We show that data on diphoton production from 
Tevatron Run II will be sensitive to a masses of the first graviton resonance in the 
range of 700-1150 GeV, while at LHC the mass range probed will be in the region of 
3.5 - 5.5 TeV. 
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The physics of extra spacetime dimensions and its imphcations for high energy physics 
and for astrophysics and cosmology has attracted tremendous attention in recent years. 
The observable universe is a dynamical hypersurface: a D^-hrane (or 3-brane) existing 
in a higher dimensional spacetime. The Standard Model (SM) fields are localized on the 
brane but gravity can propagate in the bulk. In such scenarios, the gravity /string scale 
can be lowered down from the Planck scale to the TeV scale . For high energy physics 
this is exciting because it provides fresh perspectives to the solution of the hierarchy 
problem and also suggests the discovery of new physics at TeV-scale coUiders. 

The first realization of these ideas is the ADD scenario proposed by Arkani-Hamed, 
Dimopoulos and Dvali , where, starting from a higher dimensional theory, an effective 
four- dimensional theory at a scale Ms ~ TeV is obtained. This is done by compact- 
ifying the extra dimensions to magnitudes which are large compared to the Planck 
length. The compactification radii in the ADD scenario can vary from a fermi to a 
millimetre, depending upon the number of large dimensions 0|. While the ADD model 
does not run into any obvious conflict with existing experimental data, it predicts large 
deviations from the SM in several phenomena which can be observed at existing and 
future high-energy colliders 0. Thus, laboratory data can be used to derive bounds on 
Ms in the TeV range. Bounds on this model have also been derived from cosmological 
and astrophysical considerations and some these are considerably stronger than a TeV 

i- 

The main problem that one faces within the ADD model is the reappearance of 
disparate scales viz., the string scale Ms ~ 1 TeV and the compactification radius 
Rc ~ (10~^^ TeV) ^. The stability of these Irage dimensions is an undesirable feature 
of this model and it was in an attempt to resolve this issue that the Randall-Sundrum 
(RS) model originated [^]. In its original form, the RS model is a five dimensional 
model where the fifth dimension (j) is compactified on a S^/Z^ orbifold with a radius 
Rc which is somewhat larger than the Planck length. At the orbifold fixed points, 
= 0, TT, two 3-branes called the Planck brane and the TeV brane are located. The 
SM fields are assumed to be localised on the TeV brane. To get Poincare invariance 
on the brane, it is necessary to fine-tune the cosmological constants both on the brane 
and in the bulk. The model proposes a novel five- dimensional metric of the form 

ds^ = e-'^^^^^dx^'dx" + Rld(f)\ (1) 

This metric is non-factorizable or warped and the exponential warp factor e~^^'''^ 
serves as a conformal factor for fields localised on the brane and this can be used 
to solve the hierarchy problem. The huge ratio ~ 10^^ can be generated by the 
exponent ttJCRc which needs to be only of 0(30). Rc can be stabilised against quantum 
fluctuations either by introducing an extra scalar field in the bulk H , or by invoking 
super symmetry p. 

To derive the consequences of the RS model, a linearised gravity approach is used 
where the curved metric is approximated by fluctuations h^^ about its Minkowski 
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value. On compactification of the extra dimensions, a tower of massive Kaluza-Klein 
(KK) excitations of the graviton, result on the 3-brane. The interactions of these 
with the SM particles are given by: 

iVl p iVl p ^ 

where Mp = Mp/y^Sn is the reduced Planck mass and T^"^ is the symmetric energy- 
momentum tensor for the observable particles on the 3-brane, computed using the fiat 
space metric. The masses of the h^^J are given by 

M„ = x^/Ce-"'^^^ (3) 

where the x„ are the zeros of the Bessel function Ji (x) of order unity . The masses of 
the KK excitations are not evenly spaced in this model. The zero-mode in the tower of 
excitations essentially decouples because of its weak coupling but the couplings of the 
massive RS gravitons are enhanced by the exponential e"^^^" leading to interactions 
of electroweak strength. The Feynman rules in this model are essentially the same as 



those worked out|]10], [TT| for the ADD case, except for the overall warp factor in the 
RS case. 

The basic parameters of the RS model are 

mo = /Ce""'^^^ 
Co = IC/Mp (4) 

where mo is a scale of the dimension of mass and sets the scale for the masses of the KK 
excitations, and co is an effective coupling. The interaction of massive KK gravitons 
with matter can be written as 

-T.T'''{x)h^;J{x) . (5) 



^int = - V»vr 

mo 

It is expected that the parameter cq lies in the range [0.01, 0.1]. This is because the 
scale /C is related to the curvature of the fifth dimension and so the upper bound on 
Cq results if we want to avoid strong curvature effects. But at the same time we would 
not want /C to be too small as compared to Mp, since that would introduce a new 
hierarchy. Values of mo are determined in terms of )CRc ~ 10, so that mo ranging 
from about a 100 GeV to a TeV are possible. We would like to emphasise here that 
mo cannot become arbitrarily large in the RS model. This is because if mo becomes 
very large, it would require either /C to be large, or KRc to be small (see Eq. H). This 
results in a large curvature of the fifth dimension which makes it difficult to fine-tune 
the cosmological constants on the brane and the bulk to get a fiat metric on the TeV 
brane. Consequently, the natural mass for the first graviton excitation is at most of 
the order of a few TeV. 
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It is interesting to ask what is the kind of colhder phenomenology that results 
with the RS model. Because of the fact that the zero mode decouples, it is only the 
heavier modes one can hope to detect in experiments. In the fortuitous circumstance 
that these modes are within the reach of high-energy experiments, interesting effects 
like resonance production can be observed, with the resonance decaying within the 
detectors. If this is not the case and if the the gravitons are heavier then the best 
strategy will be to look for the virtual effects of the gravitons on observables measured 
in high-energy collider experiments. Indeed, some of the phenomenology of resonant 
production of the KK excitations and the virtual effects have already been studied 
in processes like dilepton production at hadron colliders |]12[, ti production at hadron 



colliders and in deep- inelastic scattering at HERA |Tj]. Novel effects like probing 
strong gravity via black-hole production at low energies have also been discussed in 
the context of the RS model |jl5| . 

In this letter, we study the virtual effects of the exchange of spin-2 KK modes, in the 
RS model, in diphoton production at the Tevatron and the LHC In experiments, these 
photons are identified as isolated energy deposits in the electromagnetic calorimeter. 
Selection cuts are then applied to improve the signal-to-noise ratio. The DO and CDF 
collaborations have studied diphoton production in Run I of the Tevatron experiment. 
The DO collaboration has used these diphoton production data to derive constraints 
on the ADD model |jl6|. The Run I studies are however limited by statistics and 
for the purposes of the present paper of limited interest. In this paper, therefore, 
we concentrate on the production of diphotons at Run II of the Tevatron and at the 
LHC Because of the high statistics that will be achieved in these experiments, it will 
be possible to perform a detailed study of mass and angular distributions and derive 
possibly stringent constraints on the parameter space of the RS model. 

The cross-sections for the qq 77 and gg — > 77 subprocesses are O, |18[ : 



da 

-jrm 77) 

dt 



lT^C?Q\ 1 + COS^g* 

1 — cos^^* 
' ^''■Re[C(x,)](l + cos'r) + 



967r 



245767r 



(6) 



|C(x,)|^(l-cosV), (7) 



and 



da , , 
-^{99 -^11) - 

The SM box contribution gg 
energy and is negligible 



+ 6cos'r + cos^r^ 



655367r 



— s> 77 makes a very small contribution at the Tevatron 
At the LHC energy, this box contribution is somewhat 



increased because of the initial gluon flux but, as shown in Ref. [|^, in spite of this 
increase this contribution is an order of magnitude smaller than the SM qq 77 
contribution for diphoton invariant mass of 500 GeV and is more than two orders of 
magnitude smaller for diphoton invariant mass greater than about 1750 GeV. On the 
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other hand, the new physics effects dominate in the large invariant mass bins and, 
therefore, in the invariant mass region of interest the SM box contribution is negligible 
even for the case of the LHC 

In the above equations, cos^* is the scattering angle in the partonic cm. frame, 
Xs = ^ and C{x) is defined as 



Cix) = ^A(x) (9) 



with 



A(x,) = mlY^- — ^— — — . (10) 

and the Mn are the masses of the individual resonances and the r„ are the correspond- 
ing widths. The graviton widths are obtained by calculating their decays into final 
states involving SM particles. This gives 

Tn = moclxlAn (11) 

where 

A„ = Ar + Af + A^^ + A,r + ^ Ar + ^ A^ + A^^ + Af ^ (12) 

I q 

and each A^*^ is a numerical coefficient arising in the decay h"- aa. For the partial 
width A^^ , we have fixed Mh = 250 GeV in our numerical studies. 

Given the masses and the widths of the individual graviton resonances, we have to 
sum over all the resonances to get the value of X{xs)- We perform this sum numerically, 
using the fact that the higher zeros of the Bessel function become evenly-spaced. For 
a given value of = — , we retain all resonances which contribute with a significance 
greater than one per mil, and treat the remaining KK modes as virtual particles (in 
which case the sum can be done analytically). 

Using the above sub-process cross-sections we can compute the diphoton invariant 
mass distribution, da/dM and the double differential cross-section da / dMdcosO* for 
77 production cross-section at the Tevatron and the LHC, by convoluting with parton 



densities. For our numerical studies, we have used the CTEQ4M parametrisations |jT9 
for the parton distributions. To obtain the bounds on the uiq - Cq parameter space, 
we compute the cross-section for the diphoton production process in different bins of 
M (or M and cos^*, for the case of the double differential cross-section. Assuming 
Poisson-distributed data, we then compute the using the following expression: 



X^(mo,Co) = Yl 



j=bins 



„obs 

2(^f-<"^) + 2nfM-V; 



(13) 
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where n° is taken to be the SM prediction for the number of events in the bin i 
and nf^ is the prediction obtained by adding the new physics contribution to the SM 
expectation. The so calculated is used to obtain a 95% C.L. constraint on the mg - Cq 
parameter space. For Tevatron Run II, we use the following parameters: y/s = 2TeV, 

C = 2fb~\ bin size in M=80 GeV for 50 < M < 610 GeV and M = 610 1500 GeV is 

combined into one bin, y'^ < |1.2| for each of the photons. When computing the double 
differential cross-section da/dM/dcos6* , we use a bin size of 0.2 for —1 < cos^* < 1. 
For the LHC we use the following parameters: ^/s = lATeV, C = lOOfb^^, bin size 

in M=200 GeV for 500 < M < 3250 GeV and M = 3250 5000 GeV is combined 

into one bin, y'^ < |2.5| for each of the photons. NLO QCD corrections to the SM 
contribution to the diphoton cross-section yields a K-factor of 1.3 (1.1) for Tevatron 
(LHC) we assume same K-factor for the new physics contribution and simply 
multiply the full cross-section by the K-factor to account for the effect of higher-order 
QCD corrections. 




Figure 1: Constraints on the mo - Cq plane Figure 2: Constraints on the mo - cq plane 

of the Randall- Sundrum model, using 77 pro- of the Randall-Sundrum model, using 77 pro- 

duction at Run II of the Tevatron. duction at the LHC. 

We present the results of our computations in Fig. 1. In Fig 1(a), we show the 
excluded region at 95% C.L. in the (mo,Co) plane obtained by computing the as 
described above. The two curves shown in the figure are obtained by considering 
the da/dM and the da/dMdcos6* distributions. We find that the invariant mass 
distribution already gives strong constraints but the inclusion of the information on 
angular distribution does not improve these constraints. We find that for values of Cq 
between 0.01 and 0.1, the 95% C.L. limits on mo that can be derived at the Run II of the 
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Tevatron vary between 180 GeV and 290 GeV. This means that the first KK graviton 
resonance must he above 700 GeV at the least. For the LHC, the constraints are shown 
in Fig. 1(6). We find that the accessible 95% C.L. limits on mo for 0.01 < Cq < 0.1 
lie in the range between 950 GeV and 1450 GeV, which means that the first graviton 
resonance can be probed up to around 3.5 TeV. As discussed earlier, the mass of 
the first graviton excitation is at most of the order a few TeV, and it is an exciting 
prospect that this is the range of graviton masses is likely to be probed by the diphoton 
production process at the LHC. 

To summarize, we have analysed the effects of the interactions of the spin-2 Kaluza- 
Klein modes with SM fields in diphoton production at the upgraded Tevatron (Run 
II) and at the LHC, in the context of the Randall-Sundrum model. By analysing the 
invariant mass distributions of the photon pair, we derive 95% C.L. limits on rriQ for 
0.01 < Cq < 0.1. These limits are between 180 and 290 GeV for the Tevatron Run 
II and between 950 and 1450 GeV for the LHC. The range of mo values that will 
be probed by this process at the LHC are such that there is the exciting possibility 
of detecting the virtual effects at the LHC or, conversely, non-observation of these 
graviton modes at the LHC would seriously constrain the Randall-Sundrum model, at 
least in its simplest form. 
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